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Abstract. The colonization of species on remote islands may result in phenotypic diversification and ultimately speciation. 

On the Galapagos Archipelago, seven very closely related morpho-species of the wolf spiders genus Hogna are 
distinguishable based on small somatic and genital differences. Based on habitat preference, these species can broadly be 
categorized into (i) three “high elevation species” occurring on the volcanic highlands, (ii) three “coastal dry” species 
occurring in dune habitats along the coast, and (iii) one generalist species which is chiefly found in wet coastal habitats such 
as salt marshes but also in wet habitats at higher altitudes. To determine the degree of reproductive isolation among these 
morpho-species, we investigated gene flow among populations and species based on nine allozyme loci. Genetic analysis by 
means of genetic distance estimates and cluster agglomerative analyses confirmed the status of the defined morpho-species. 

Allele frequencies were highly similar among populations within a species but differed profoundly among species. Genetic 
differentiation within the generalist speeies was generally very low. There were no constant differences between high 
elevation and coastal populations for this species. Neutral genetic divergence between species appeared to correspond more 
to geographic distribution rather than to a clear separation of the two different ecological groups within an island. This 
suggests that a parallel parapatric divergence between high elevation and coastal dry species may have taken plaee on the 
oldest islands of San Cristobal and Santa Cruz. 

Keywords: Wolf spiders, speciation, island biogeography, ecological speciation 


Archipelagos are among the world’s great natural labora¬ 
tories of evolution, as many studies on the Galapagos, 
Hawaiian, Canary Islands, and other island groups have 
shown. The Galapagos are of particular interest for the 
following reasons: they are truly volcanic, well isolated 
(between 900 and 1000 km west of the Ecuadorian mainland), 
and of known age (Simkin 1984). There is no evidence of the 
existence of land bridges so all terrestrial organisms had to 
cross an oceanic barrier by dispersal from the mainland. 

The Galapagos Archipelago consists of 13 large islands and 
a great number of islets and rocks, all of volcanic origin 
(Fig. 1). The southeastern islands are the oldest (3-5 million 
years) while the northern and western islands are the youngest 
(< 0.7 million years) (Simkin 1984). 

Due to geographic isolation, many endemic animal (e.g., 
Darwin’s finches, giant tortoises, lava lizards, mockingbirds) 
and plant (e.g., Opuntia cacti, Sccilesia trees) groups have 
radiated. Evolutionary research on these islands has mainly 
focused on vertebrate species such as Darwin’s finches, giant 
tortoises, lava lizards, mockingbirds, and on plant species such 
as Opuntia cacti and Scalesia trees (Grant 1981; Fritts 1984; 
Snell et al. 1984; Stern & Grant 1996; Rassmann et al. 1997; 
Cacone et al. 2002). Speciation patterns of invertebrates have 
been, in contrast, much less studied, and only recently have 
genetic studies been conducted on Coleoptera such as the 
tenebrionid Stomium (Finston Sc Peck 1995), the chrysomelid 
Nesaecrepidia (Verdyck Sc Desender 1999), the carabid Calo- 
soma (Desender & Verdyck 2000; Verdyck et al. 2003,2004), the 
weevil Galapaganus (Sequeira et al. 2000, 2008), and the land 


snail genus Bulimnlns (Parent & Crespi, 2006). Genetic studies 
on Galapagos spiders are presently lacking, while such studies 
in other locations have revealed adaptive radiations on other 
archipelagos (e.g., Hawaian Tetragnatha spiders, Gillespie 
2004; Hawaiian Dysdera spiders, Arnedo 2001). 

Previous studies on the spider genus Hogna , the only wolf 
spider genus occurring on the archipelago (Maelfait & Baert 
1986; Baert & Maelfait 1997), revealed that this genus consists 
of several closely related, or even cryptic, species. Based on 
somatic and small genital differences, a total of seven morpho- 
species is suggested with a distinct distribution on the islands 
(Baert et al. unpublished data) (Fig. 1). At least three groups 
of morpho-species can be distinguished that differ distinctly. 

A first group of ecologically and morphologically similar 
morpho-species occur at higher altitudes on the islands in the 
pampa vegetation zone and are hereafter referred to as “high 
elevation species.” Based on differences in morphology of the 
genital organs, different species can be distinguished: Hogna 
species 1 (HI), living on both southern volcanoes Cerro Azul and 
Sierra Negra of Isabela, one of the youngest islands, Hogna 
species 4 (H4) which occurs on islands of intermediate age (Santa 
Cruz, Santiago, and Volcan Alcedo of Isabela), and Hogna 
species 2 (H2) whieh occurs on the oldest island of San Cristobal. 

A second group, hereafter referred to as “coastal dry 
species,” lives in the dry arid zone along the coast in vegetated 
dunes and in the Opuntia cactus zone. These morpho-species 
can only be found on the oldest islands of San Cristobal 
(Hogna species 5 (H5)), Espanola (Hogna species 7 (H7)), and 
Santa Cruz (Hogna species 6 (H6)). The San Cristobal species 
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Figure 1.—Map of the sampling locations of the Hogna species from Galapagos. Nearby sampling locations are presented as one dot due to 
space restrictions. Dots between brackets for H4 indicate sites where the species is present but that were not included in this study. 


(H5) is mostly found in the depressions overgrown with sea 
grass (Sporobohis virgimcus) located behind the shore or on 
low vegetated dunes, while H7 is found in tall vegetation of the 
dry arid zone directly adjacent to the littoral zone, but never in 
the adjacent depressions with salt grass. The Santa Cruz 
species (H6) is found in the Opimtia cactus zone in between 
dune and pure rocky soil. 

The third group comprises populations of the generalist 
species Hogna species 3 that lives in saline habitats along the 
coast (salt marshes, bays), along permanent pools (e.g., El 
Chato on Isla Santa Cruz) and in permanent wetlands below 
600 m of altitude (e.g., Los Gemelos on Isla Santa Cruz). 
Scattered populations can also be found above the vegetation 
inversion zone in wet conditions during El Nino years 
(characterized by very heavy rainfall giving rise to temporary 
pools) (Baert & Maelfait 2000). They reach, however, their 
highest densities in the salt marshes. All these populations 
have very similar genital organs and are at present interpreted 
as belonging to a single species. It is widespread over the whole 
archipelago, with the exception of the northern island Pinta 
and the southeastern island Espanola (Baert & Maelfait 1997). 

In this paper, we test whether the separation of this genus 
into seven morpho-species on the Galapagos is justified. By 
means of cellulose acetate gel electrophoresis, wherein 8 
allozyme loci (FUM, G6PDH" GOT, IDH, LDH, MPI, PGI, 
and PGM) were studied, we investigate whether the genetic 
variation among species is larger compared to the variation 
among populations within species and indicative of reproduc¬ 
tive isolation among the species. 

METHODS 

Sampling collection*— In the period between 1996 and 2002, 
we sampled a total of 43 known Hogna populations (see 
Table 1) from 9 islands (Santa Cruz, Isabela Volcan Sierra 


Negra, Isabela Volcan Cerro Azul, San Cristobal, Floreana, 
Rabida, Genovesa, Bartolome, Santiago, and Espanola) and 
seven morpho-species. In three localities, the high elevation 
species occurred sympatrically with H3, [e.g., Cerro Gavillan 
(populations 40 & 41), El Junco (populations 1 & 2), and Los 
Gemelos (population 23) (Fig. 1)]. Populations of H4 occur¬ 
ring on the tops of Isabela and Santiago were preserved in 
ethanol and could therefore not be included in this allozyme 
study. 

Individuals were caught by hand, mostly at sunset with an 
electric torch worn on the forehead. They were stored and 
transported in a Taylor-Wharton cryogenic shipper saturated 
with liquid nitrogen. In the laboratory, the material was stored 
in an ultra-cold freezer at —80° C. The aim was to investigate 
at least 40 individuals for each population if possible. In some 
localities, their densities were so low that this number could 
not be reached. Some localities were sampled several times but 
in different years. Voucher specimens are deposited at the 
Royal Belgian Institute of Natural Sciences. 

Allozynics.—Parts of the legs were homogenized in distilled 
water for performing the cellulose acetate gel electrophoresis, 
following the procedures of Hebert & Beaton (1989). Eight 
enzymes (9 loci) were tested for polymorphism: fumarate 
hydratase (FUM), aspartate minotransferase (GOT1, GOT2), 
isocitrate dehydrogenase (IDH),lactate dehydrogenase 
(LDH), mannose phosphate isomerase (MPI), phosphoglu- 
cose isomerase (PGI), phosphoglucomutase (PGM), and 6- 
phosphogluconate dehydrogenase (6PGDH). 

Allele frequencies were obtained for each population and 
species. Deviations from Hardy-Weinberg equilibrium were 
tested by means of an exact test. Genetic divergence between 
populations and species were estimated based on Nei’s 
unbiased genetic distance (Nei 1978). Based on this distance 
metric, divergence among populations within species was 
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Table 1.—List of the 43 sampled Hogna populations on Galapagos during the years 1996, 1998, 2000, and 2002. Situation and number of 
caught specimens per sample. Abbreviations: SCB = Isla San Cristobal; ESP = Espahola; GEN = Genovesa; SCZ = Santa Cruz; FLO = 
Floreana; BAR = Bartolome; RAB = Rabida; SAN = Santiago; ISN = Isabela, Volcan Sierra Negra; ICA = Isabela, Volcan Cerro Azul. 


Code 

Island 

Locality 

Vegetation zone 

Elevation 

Morpho-species 1996 

1998 

2000 

2002 

1 

SCB 

El junco 

pampa 

675m 

H3 

8 



36 

2 

SCB 

El junco 

pampa 

675m 

H2 

3 



44 

3 

SCB 

Cerro San Joaquin 

pampa 

700m 

H2 




46 

4 

SCB 

Punto Baso 

dune 

5 m 

H5 




41 

5 

SCB 

Punto Baso (fregat nesting) 

dune 

5 m 

H5 




36 

6 

SCB 

Punto Baso (Sesuvinm) 

littoral zone 

lm 

H3 




27 

7 

SCB 

Caleta de la Tortuga 

littoral zone 

5m 

H3 




10 

8 

SCB 

La Loberia 

littoral zone 

lm 

H3 




42 

9 

SCB 

Caleta Sapho (Spirobolus) 

salt grass 

lm 

H5 

28 



55 

10 

ESP 

Punta Cevallos 

dry arid zone 

2m 

H7 



40 


11 

ESP 

Bahia Gardner 

dry arid zone 

lm 

H7 



31 


12 

ESP 

Isla Gardner 

dry arid zone 

2m 

117 



41 


13 

GEN 

Lago Arcturus 

littoral zone (lagoon) 

60 m 

H3 



18 


14 

SCZ 

Laguna Andreas 

littoral zone (lagoon) 

lm 

H3 

35 


42 

40 

15 

SCZ 

Bowdich 

littoral zone (lagoon) 


H3 




31 

16 

SCZ 

El Garapatero 

littoral zone (lagoon) 


H3 



58 

22 

17 

SCZ 

Las Palmas 

dry arid zone 


H6 




28 

18 

SCZ 

Meteo Station 

littoral zone (lagoon) 

5m 

H3 


40 



19 

SCZ 

Bahia Tortuga 

littoral zone (lagoon) 

lm 

H3 


41 

27 

53 

20 

SCZ 

Playa Bachas 

littoral zone (lagoon) 

lm 

H3 


43 

40 

44 

21 

SCZ 

El Chato 

around permanent pool 


H3 



40 

40 

22 

SCZ 

El Carmen 

temporary pool (El Nino) 

350m 

H3 


42 



23 

SCZ 

Los Gemelos, open 

pampa 

570m 

H3 




43 

24 

SCZ 

Los Gemelos, Scalesia 

Scalesia forest 


H3 



47 

40 

25 

SCZ 

Media Luna 

pampa 

600m 

H3 


51 

37 

27 

26 

SCZ 

Tss ML & Cpunt 

pampa 


m 




13 

27 

SCZ 

Cerro Crocker 

pampa 

875m 

H4 


42 

23 

47 

28 

FLO 

Punta Cormoran 

littoral zone (lagoon) 

lm 

H3 



40 


29 

FLO 

Finca Cruz 

pampa 

200m 

H3 



15 


30 

FLO 

Highland 

pampa 

350m 

H3 



50 


31 

BAR 


littoral zone (lagoon) 

lm 

H3 


31 



32 

RAB 


littoral zone (lagoon) 

lm 

H3 


38 



33 

SAN 

Playa Espumila 

littoral zone (lagoon) 

lm 

H3 


45 


40 

34 

SAN 

Aguacate 

transition zone (El Nino) 

500m 

H3 


40 


41 

35 

SAN 

La Central 

pampa (El Nino) 

700m 

H3 


38 


40 

36 

SAN 

Jaboncillo 

pampa (El Nino) 

820m 

H3 


40 


40 

37 

ISN 

Laguna de Villamil 

littoral zone (lagoon) 

lm 

H3 

41 




38 

ISN 

Top 

pampa 


HI 

42 




39 

ICA 

Caleta Iguana 

littoral zone/dry arid zone 

2m 

H3 


32 



40 

ICA 

Cerro Gavilan 

pampa 

700m 

H3 


12 



41 

ICA 

Cerro Gavilan 

pampa 

700m 

HI 


30 



42 

ICA 

1100m 

dry arid zone (El Nino) 

1100m 

H3 


11 



43 

ICA 

Top 

dry arid zone (El Nino) 

1530m 

H3 


30 





Annual total no. specimens 




157 

606 

549 

926 


compared with the distance among populations of different 
species. These analyses were performed with the computer 
packages TFPGA (Miller 1997) and GenAlEx (Peakall and 
Smouse 2006). Genetic distances were visualized by means of 
principal component analysis (PCA), designed for ordination 
of allelic frequency data, by means of the computer package 
PCA-Gen (Goudet 1999). 

RESULTS 

Allelic variation and heterozygosity were very low within 
each species, but differed clearly among species, with one or 
a few alleles that were fixed within the morpho-species. The 
low genetic variability among populations within species, 
compared to the variability among species, is clearly 


depicted when genetic distances are compared among 
populations (Table 2). The genetic distance between popu¬ 
lations belonging to the same morpho-species ranged from 0 
to maximum 0.031 (H3), demonstrating that the allele 
frequencies of the different populations within a given 
morpho-species were highly similar. Differences in allele 
frequency of populations belonging to a different morpho- 
species were, in contrast, considerably higher and ranged 
from 0.118 to 2.277. The smallest genetic distances were 
observed between H2 (“high elevation” San Cristobal) and 
H7 (“coastal dry” Espahola) and between H5 (’’coastal dry” 
San Cristobal) and H7. 

Allele frequencies for all loci were near to fixation for 
almost all morph-species (Table 3). None of the species 
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Table 2.—Genetic distance data for seven Hogna morpho-species on the Galapagos. Above diagonal: genetic distances between the different 
morpho-species when allele frequencies of the different populations of each morpho-species were pooled; diagonal: minimum and maximum 
genetic distance between populations of the same morpho-species; below diagonal; minimum and maximum genetic distance between 
populations of different morpho-species. 


HI 

H2 

H3 

H4 

H5 

H6 

H7 

HI 0.012 

1.711 


1.026 

0.449 

1.218 

1.089 

1.719 

H2 1.647-1.811 

0.0016 


2.059 

0.587 

0.267 

1.151 

0.131 

H3 0.910-1.040 

1.614-2.154 

0 . 000 - 

-0.031 

1.481 

1.474 

1.090 

2.193 

H4 0.435-0.483 

0.577-0.600 

1.230 

1.504 

0.0001 

0.577 

0.603 

0.800 

H5 1.179-1.281 

0.260-0.275 

1.213- 

1.503 

0.575-0.584 

0-0.0001 

0.817 

0.125 

H6 1.084-1.097 

1.123-1.182 

1.083- 

-1.235 

0.602-0.607 

0.813-0.822 

— 

1.119 

H7 1.645-1.841 

0.122-0.146 

1.727- 

-2.277 

0.798-0.811 

0.118-0.138 1.121-0.122 

0.000-0.013 

comparisons resulted in fixation of the same alleles, and each from PCA ordination of the different populations (Fig. 2). 

morpho-species was, therefore 

, characterized by a unique The first three axes explained 74.24%, 15.06%, 

and 4.06% of 

allele combination. 



the total allelic variation respectively. Along the first axis, 

The low differences among con specific populations 

com- the H3 populations are clearly separated from the other 

pared to differences among morpho-species were also obtained morpho-species. The position of the remaining species along 

Table 

3.- Allele frequency data for seven Hogna morpho-species on the Galapagos. 




HI 

H2 

H3 

H4 

H5 

H6 

H7 

Npop 

2 

2 

30 

2 

3 

1 

3 

Nind 

72 

91 

1632 

144 

112 

28 

153 

FUM 1 

0.0000 

0.0000 

0.0000 

0.0035 

0.0000 

1.000 

0.0000 

2 

0.9931 

0.0000 

0.0000 

0.9896 

0.0000 

0.0000 

0.0000 

3 

0.0000 

0.9890 

0.0000 

0.0035 

1.0000 

0.0000 

1.0000 

4 

0.0069 

0.0110 

1.0000 

0.0035 

0.0000 

0.0000 

0.0000 

GOT1 1 

0.0000 

0.0165 

0.0000 

0.0000 

0.0000 

0.0000 

0.0033 

2 

0.0278 

0.9560 

0.0018 

1.0000 

1.0000 

1.0000 

0.9967 

3 

0.9722 

0.0220 

0.9982 

0.0000 

0.0000 

0.0000 

0.0000 

4 

0.0000 

0.0055 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

LDH 1 

1.0000 

0.0000 

0.9997 

0.0035 

0.0000 

0.0000 

0.0000 

2 

0.0000 

0.0000 

0.0000 

0.0069 

1.0000 

0.0000 

1.0000 

3 

0.0000 

1.0000 

0.0003 

0.9861 

0.0000 

0.0000 

0.0000 

4 

0.0000 

0.0000 

0.0000 

0.0035 

0.0000 

1.0000 

0.0000 

G6P3H 1 

0.0000 

1.0000 

0.0003 

0.0000 

0.0000 

0.0000 

1.0000 

2 

1.0000 

0.0000 

0.9994 

0.9931 

1.0000 

1.0000 

0.0000 

3 

0.0000 

0.0000 

0.0003 

0.0069 

0.0000 

0.0000 

0.0000 

MPI 1 

0.0000 

0.0714 

0.9997 

0.0035 

0.0045 

0.0000 

0.0000 

2 

0.0000 

0.9286 

0.0003 

0.0070 

0.9955 

0.0000 

1.0000 

3 

1.0000 

0.0000 

0.0000 

0.9860 

0.0000 

0.0000 

0.0000 

4 

0.0000 

0.0000 

0.0000 

0.0035 

0.0000 

1.0000 

0.0000 

PG1 1 

0.0000 

0.0330 

0.0003 

0.0139 

0.0000 

0.0179 

0.0000 

2 

0.0069 

0.8791 

0.9988 

0.9549 

0.9911 

0.9286 

0.9412 

3 

0.0139 

0.0824 

0.0009 

0.0174 

0.0089 

0.0536 

0.0588 

4 

0.9583 

0.0000 

0.0000 

0.0139 

0.0000 

0.0000 

0.0000 

5 

0.0208 

0.0055 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

PGM 1 

0.0000 

0.0275 

0.9969 

0.0035 

0.0000 

0.0536 

0.0000 

2 

0.0278 

0.9341 

0.0031 

0.0000 

0.9688 

0.0000 

0.8170 

3 

0.9722 

0.0330 

0.0000 

0.9931 

0.0313 

0.9464 

0.0065 

4 

0.0000 

0.0055 

0.0000 

0.0035 

0.0000 

0.0000 

0.1765 

GOT2 1 

0.0625 

0.0495 

0.9991 

0.0069 

0.0089 

0.0000 

0.0000 

2 

0.9375 

0.9286 

0.0009 

0.9861 

0.9911 

0.9286 

1.0000 

3 

0.0000 

0.0220 

0.0000 

0.0069 

0.0000 

0.0714 

0.0000 

IDH1 1 

0.0000 

0.0000 

0.9694 

0.0104 

0.0000 

1.0000 

0.0000 

2 

0.5000 

1.0000 

0.0306 

0.9896 

1.0000 

0.0000 

1.0000 

3 

0.5000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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Figure 2.—Results of principal component analysis (axes 1, 2 and 3). (Filled squares = high elevation species; filled diamonds = coastal dry 
species; open circles = generalist species.) 


this axis and the second axis corresponded to their 
geographic position and age of the islands rather than to 
their habitat preference. The three species of the oldest 
islands - H7 from Espanola and H2 and H5 from San 
Cristobal - were clustered in the PCA, followed by H4 and 
H6 from Santa Cruz and HI from the volcanoes Cerro Azul 
and Sierra Negra on Isabela. 

DISCUSSION 

Our results show that there is a clear and very high degree of 
genetic divergence between the previously defined morpho- 
species. Moreover, genetic divergence among populations 
within species was much lower (Table 2, Fig. 2). These two 
findings indicate that these morpho-species likely represent 
distinct reproductively isolated species. 


Although the validity of allozymes for phylogenetic 
inferences is questionable since the historical and genealogical 
relationship between the different alleles remains unknown 
(Lowe et al. 2004), a few suggestions concerning historical 
patterns of divergence can be deduced. 

First, these results suggest that Hogna speciation on the 
Galapagos is likely due the combined effect of geographic 
isolation and ecological specialization within the different 
climatological and vegetation zones present on the different 
islands. Except for H3, species from the same or proximate 
island tend to be genetically more similar to each other. 
Species living on the same islands but in different habitats are, 
however, genetically fixed for a few loci, clearly indicating a 
lack of gene flow and hence strong reproductive isolation. 
Combining these results suggests that ecological specialization 
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on the islands Santa Cruz and San Cristobal occurred 
repeatedly in association with speciation events rather than a 
diversification of a habitat adapted lineage with secondary 
colonization of the specialized forms to the different islands. 
According to the second scenario, species living in the same 
habitat would then be expected to be genetically more similar 
to each other. 

Similar patterns of species diversification in terms of 
geographic position and ecological specialization have been 
confirmed by more thorough genetic analyses of Tetragnatha 
spiders from Hawai (Gillespie 2004). In the Galapagos, this 
speciation pattern has been observed in other terrestrial 
invertebrates such as weevils and snails (Parent & Crespi 
2006; Sequiera et al. 2008). 

Whether the generalist species H3 can be regarded as closely 
related to the ancestral species, as suggested by Baert & 
Maelfait (2000), however, cannot be confirmed by these data. 
The smaller genetic distance between the generalist 113 and the 
species from younger islands (HI) compared to those of older 
islands (H2, H5, and H7) is in accordance with this hypothesis. 
Surprisingly however, H3 showed a very low degree of genetic 
variation within populations. Moreover, distant populations 
as well as populations living in different habitats all appeared 
to be genetically very homogenous. These observations 
suggest that this generalist and apparently highly dispersive 
species may have colonized the archipelago independently. 

The results can only be interpreted as preliminary as they 
are based on allozyme data and only a few loci were scored. 
Moreover, the selective neutrality of allozymes has often been 
questioned. Our ongoing work aims to add more variable loci 
such as mitochondrial DNA (Cytochrome Oxidase I) so that 
more well founded phylogenetic inferences can be made. Also, 
future work should include Hogna species from the South 
American mainland to better understand the phylogenetic 
relationships between the species and the colonization history 
of Hogna in Galapagos. 
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